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The mass spectrum of heavy quark baryons has been computed on anisotropic lattices using quenched lattice 
nonrelativistic QCD. The mass splittings between spin-1/2 and spin-3/2 baryons are also calculated. Results are 
compared to those obtained by using a Dirac- Wilson action of the D234 type. Color hyperfine effects in heavy 
baryons are also discussed. 



1. INTRODUCTION 

Extraction of correct mass splittings between 
vector and pseudoscalar mesons remains a chal- 
lenging problem in lattice QCD. While quenched 
lattice QCD has been unable to extract correct 
mass splittings so far this problem persists 
even with full lattice QCD |^ . One may thus ask 
whether lattice results for baryon mass splittings 
exhibit similar suppression compared to experi- 
mental numbers. So far, only a few lattice QCD 
results have been reported for heavy baryons 
, and there are only two works where heavy 
quarks are treated nonrelativistically. In this 
work we report heavy baryon results by using 
a nonrelativistic heavy quark action and an im- 
proved light quark action on anisotropic lattices. 

2. NUMERICAL SIMULATION 

2.1. Actions 

For light quarks we use a Dirac- Wilson action 
of the D234 type ||], which has been used pre- 
viously and detailed in Refs. H,^. Its leading 
classical errors are cubic in lattice spacing. The 
gauge action as well as the heavy quark NRQCD 
action is detailed in Ref. The gauge action 
is tadpole improved and the leading classical er- 
ror is quartic in lattice spacing. The Hamiltonian 
corresponding to the NRQCD action is complete 
to 0{1/M^) in the classical continuum limit. 

2.2. Simulation Details 

This work is done with two sets of quenched 
gauge configurations on anisotropic lattices (lat- 



tice sizes 12^ X 32 and 14^ x 38 at /3 = 2.1 and 2.3 
respectively) with a bare aspect ratio as/ at = 2. 
For (3 ~ 2.1, we use 720 configurations and for 
(3 ~ 2.3, the number of configurations is 442. 
Two sets of bare masses are used for heavy quarks 
while four sets of hopping parameters are used for 
light quarks and they are listed in Table 1. In- 
terpolating operators used here are same as in 
Ref. [^,0, and correlation functions are calcu- 
lated using both local and smeared (gauge invari- 
ant smearing at the sink ^,^,0) propagators. The 
time window to extract mass is chosen in a way 
such that the ending time is large and the fit is 
stable with variation of both starting and ending 
time by a few time steps. The fit corresponding to 
lowest is considered in this window-span. The 
temporal lattice spacing and correspondingly the 
scale is fixed by extracting the p-meson mass and 
the strange mass is fixed from the (/)-meson mass 
. The charm mass is fixed by extracting the rj^ 
mass, whereas B'^ mass is used to fix the bottom 
mass. As in Ref. light quark extrapolation 
is done by extrapolating the hadron masses as a 
function of cq -I- C2»n^ -l- csto^, where is the 
pion mass. 

The kinetic mass of a state can be extracted by 
using the usual NRQCD relation 

27r2 

^'^^ = m,eat[at{E,-E,)y 

which is derived from E = ^^^^ (see Ref. |^] 
for details and clarification of notations). Mass 
differences between two states [Hi and H2) with 
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Table 1 

Summary of lattice parameters. 



/3 



K (light quark) 



bare mass 
c h 



2.1 0.229,0.233,0.237,0.240 1.2,1.5 5.0,6.0 
2.3 0.229,0.233,0.237,0.240 1.04,1.24 3.7,4.2 



the same heavy quark can be obtained by taking 
the difference of their zero momentum simulation 



energy : 



(2) 



It is always easier (smaller error) to extract heavy 
meson or baryon masses with heaviest light quark 
{k = 0.229) rather than extracting their masses 
with lighter light quark {k = 0.233 and higher). 
Moreover, we can calculate mass differences more 
accurately rather than masses, and so one can use 
mass differences to a extract mass. For example, 
one can calculate the meson mass by using the 
relation 



M{qi,Q) = Miqh,H)-AM 
where 

AM = AE ^ E{qh,Q) - Eiqi,Q). 



(3) 



(4) 



Here qi and qh denote lighter and heaviest light 
quark respectively. This equation is valid as long 
as the heavy quark Q is same in both states. Sim- 
ilarly, masses of singly and doubly heavy baryons 
can be extracted from meson masses by using 



M(<Zi92,Q) = 
M{qi,QQ) = 

where 



M{qh,Q)-AE,h, 
M{QQ) - AEdh, 



Esh = E{qh,Q)~E{qiq2,Q), 
Edh = E{QQ)-E{qi,QQ). 



(5) 
(6) 

(7) 

(8) 



For example, Sc(b) mass is extracted by taking 
its difference (at each k) with D{B^) mass (m) 
a.t K = 0.229 and then subtracting that from m. 
Masses extracted by using Eq. 1 and Eqs. 3-8 are 
consistent with each other. However, errors in the 
second method are considerably smaller than the 
previous one. 



Table 2 

NRQCD Results (in units of MeV). First error 
is statistical while second error comprises sys- 
tematic errors due to scale, time window and 
anisotropy. 



(3 = 2.1 



'CC 



2407(32)(|?)" 
2652(25)(2I) 
75(20)(i|) 
71(18)(P) 
65(13) (D 
128(28)(i) 

"3562(47)(ir 
3681(44)(}D 
63(14)(?) 
56(8)(D 



2.3 

2471(38X1)" 
2690(33) (ti) 
90(18)(}|) 
86(16)(}i) 

78(14)(?o) 
155(38)(i) 



- ^b 
"b ~ "fc 

ni — rib 

Ab-Sb 



^bb ^bb 
^bb ~ ^bb 
^cb 

' 'cb ' 'cb 

O* — O' 

"cb "cb 
^cb — ^cb 



cb 



cb 



22(12)(i) 
21(12)(i) 
18(8)(!) 
142(28)(i) 

22(6)(^) 
20(4)(i) 
6810(150)(?^) 
6935(135)(i) 
46(8)(|) 
40(6)(|) 
11(6)(|) 
10(5)(^) 



3580(67)(|J) 
3692(62)(p) 
72(11)(D 
65(7)(I) 
5681(102)(tg) 
26(13)(i) 
26(12)(i) 
24(8)(i) 



21(6)(^) 
20(4)(i) 
6817(230)(f^i 
6928(217)(' 
48(9)(^ 
43(6)(i) 
10(5)(1) 
10(4)(1) 



('62 , 
115^ 



3. RESULTS AND SUMMARY 

The mass spectrum and spin splittings of 
heavy quark baryons have been computed on 
anisotropic lattices using an NRQCD heavy quark 
action. Results are summarized in Table 2, 
whereas in Table 3 we have compared our results 
with those obtained by using a relativistic heavy 
quark (D234 action, ^) and experimental num- 
bers (where available). One can notice NRQCD 
results and D234 results are consistent with each 
other. Results are also consistent with a previous 
NRQCD calculation @. 

As it has been argued earlier p,0, suppres- 
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Table 3 

Comparison of NRQCD results with D234 re- 
sults and experiment [0. For lattice results, 
first row is for D234 while the second row is for 
NRQCD results (all are in MeV). 





Lattice 


results 


Expt 


0.4 








0.3 




f3 = 2.1 


(3 = 2.3 






2379(31)(^^) 


2490(14)(^^) 


2455 






2407(32)(P) 


2471(38)(i9) 




0.2 




2671(ll)(ii) 


2699(10)(|i) 


2704 


0.1 




62(33)(i) 
75(20)(1|) 


82(12)(i) 
90(18)(}|) 


64 






52(15)(1) 
71(18)(i2) 


82(10)(i) 
86(16)(}i) 


70 


0.4 


n* - nc 


50(17)(ii) 


73(8)(D 




0.3 




65(13)a) 


78(14)(?o) 






^cc 


3608(15)(Jf) 
3562(47)(i) 


3595(12)(2i) 
3580(67)(|^) 




0-2 


^cc 


3747(9) (1^) 
3681(44)(}^) 


3727(9)(16) 
3692(62)(i^) 




0.1 




58(14)(18) 


83(8)(7o) 









63(14)(?) 


72(11)(D 
72(5)(t) 








57(8)(i°) 








56(8)(D 


65(7)(I) 
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splittings as a function of average meson mass. 



sion of spin splittings is not present in baryon 
sector in the same way as it is in heavy meson 
sector (characteristic of quenched QCD simula- 
tions). In Fig. 1 , we plot meson and baryon 
spin splittings respectively. There is a clear sup- 
pression in meson sector (top fig) which is not 
present in baryon sector (bottom fig). Exper- 
imental measurement of spin splittings in bot- 
tom baryons as well as for doubly heavy baryons 
will be very useful to check if this is true indeed. 



This work is supported by the Natural Sciences 
and Engineering Research Council of Canada. 

REFERENCES 

1. For a review see C. Bernard, Nucl. Phys. 
B(Proc. Suppl.) 94, 159(2001). 

2. C. Stewart and R. Koniuk, Phys. Rev. D 63, 
054503(2001); T. Manke et al, Phys. Rev. D 



62, 114508(2000); S. Collins et at, Phys. Rev. 
D 60, 074504(1999). 

3. UKQCD Collaboration, K.C. Bowler et al, 
Phys. Rev. D 54, 3619(1996). 

4. A. Ah Khan et al., Phys. Rev. D 62, 
054505(2000). 

5. R.M. Woloshyn, Phys. Lett. B 476, 
309(2000). 

6. Randy Lewis, Nilmani Mathur and R.M. 
Woloshyn, hep-ph/0107037. 

7. Randy Lewis, Nilmani Mathur and R.M. 
Woloshyn, hep-lat/0109014. 

8. M. Alford, T.R. Klassen, and CP. Lepage, 
Nucl.Phys. B496, 377(1997). 

9. R. Lewis and R.M. Woloshyn, Phys. Rev. D 
62, 114507(2000). 

10. D.E. Groom et al, Eur. Phys. J. C 15, 
1(2000); CLEO Collaboration, CP. Jessop et 
al, Phys. Rev. Lett. 82, 492(1999). 



